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NOMBRE DE RENCONTRES ENTRE LES INDIVIDUS 
DE N POPULATIONS ANIMALES HOMOGENES. 
ETABLISSEMENT D'UNE EQUATION. 

PAR 

R. R. Tercafs ( 1 ). 


INTRODUCTION. 

Le but général de notre travail est la détermination des facteurs régis¬ 
sant la formation, l'évolution et le maintien des populations animales. 
Ce programme comprend, d'une part, la transformation de certaines don¬ 
nées biologiques essentielles en paramètres mathématiques, d'autre part, 
l'élaboration de modèles mathématiques exprimant les relations entre les 
divers composants d'une population. 

L'évolution d’une population animale est la résultante d'un grand 
nombre de facteurs. Cependant, certains d'entre eux sont prédominants. 
Parmi ceux-ci, on peut citer : le nombre de proies détruites par les car¬ 
nivores, le taux de reproduction des différents biotes, les masses méta¬ 
bolisées, etc. 

Dans la plupart de ces phénomènes apparaît une donnée essentielle : 
le nombre de rencontres entre les individus. Dans le cas des prédateurs, 
il est bien évident que le nombre de proies détruites, par unité de temps, 
dépendra du nombre de contacts entre les animaux appartenant à ces 
deux niveaux trophiques. De même, le taux de reproduction sera, égale¬ 
ment, en partie fonction du nombre d'accouplements, c'est-à-dire du 
nombre de rencontres entre les individus de sexe différent d’une même 
population. 

Le nombre de rencontres entre les individus d'une population est donc 
un processus essentiel régissant l'évolution de cette population. 


(I) Chargé de recherches du Fonds National de la Recherche Scientifique de Bel¬ 
gique (F. N. R. S.) 

Bénéficiaire d'un crédit aux chercheurs du F. N. R. S. 





2 


R. R. TERCAFS. — RENCONTRES ENTRE INDIVIDUS 


44, 36 


Toute élaboration de modèle mathématique tendant à exprimer de 
façon précise les variations des divers constituants d'une population ani¬ 
male suppose donc la quantification préalable de cette relation fonda¬ 
mentale. 

Nous avons mis au point une équation exprimant le nombre de rencon¬ 
tres entre les individus de n catégories différentes en fonction de divers 
paramètres. C’est le détail de cette formulation, avec des exemples d’expé¬ 
riences, que nous présentons dans ce travail. 


Etablissement d’une équation 
exprimant le nombre de rencontres 
entre les individus de n populations différentes. 

L’équation sera d’abord établie pour deux populations puis étendue 
à n populations. 

Considérons une espace E contenant une population animale PI, distri¬ 
buée dans tout l'espace de façon aléatoire. Supposons que les animaux de 
cette population se déplacent tous, pendant l'espace de temps considéré. 

Cette population est caractérisée par les facteurs suivants: 

m = nombre d'animaux de la population Pl par cm 2 . 

Dl = diamètre d’influence ou de perception d'un animal (2). 

VI = vitesse moyenne de déplacement des animaux de Pl. 

Considérons, d'autre part, un animal, A(P2), appartenant à une autre 
population P2, vivant dans le même espace E. Un animal de la popu¬ 
lation P2 est caractérisé par : 

V2 = vitesse moyenne de déplacement. 

D2 = diamètre d’influence ou de perception. 

Dans un cm 2 de l’espace E, il y a m animaux Pl. Les animaux se 
déplacent dans le même plan. Il y aura rencontre entre l’animal A(P2) 
et un animal de Pl si les deux surfaces d’influence ou de perception 

TT D2 2 ttDI 2 

- et - viennent en contact. Ou plus simplement, puisque le 

2 2 

mouvement a lieu dans un plan, lorsque les diamètres Dl et D2 se ren¬ 
contrent. La longueur totale de contact possible entre l’animal 
A (P2) et les animaux Pl est donc égale à m D, avec D = Dl + D2. 

(2) Dans le cas des objets inertes, on ne fera intervenir que le diamètre réel de l'objet. 
Dans le cas des animaux, il est bien évident qu'il n'y a pas besoin d'un contact physique 
pour qu'il y ait rencontre entre deux individus. On doit alors faire intervenir soit un 
diamètre de perception (dans le cas d’un prédateur, par exemple) ou d'influence : une 
proie émet autour d'elle un certain nombre d'informations (odeur, substances sur le sol, 
vibrations, etc.). 
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D'autre part, chaque animal de la population PI se déplace dans 
toutes les directions avec une vitesse VI. L'animal A(P2) se déplace 
dans toutes les directions avec une vitesse V2. Il existe donc une vitesse 
de déplacement relative entre deux animaux de PI et de P2 pris 
deux à deux. Appelions cette vitesse Vr. 

En une seconde, l’animal A(P2) s'est déplacé par rapport aux animaux 
Pl d'une distance dont la valeur est donnée par Vr. 

La longueur totale de contact possible est donc maintenant égale à 
m D Vr, Vr étant la vitesse relative moyenne. 

Le nombre de rencontres, par seconde et par cm-, entre u n animal de 
P2 et les animaux de la population Pl est donc: 

= m D Vr (1). 

S’il y a n animaux de P2, on a : 

0 12 = n m D Vr (2). 

De même, dans un système à n populations, on calculera le nombre 
de rencontres pour chaque catégorie prises deux à deux. 

Calcul de Vr. 

Les directions de déplacement de deux animaux A(P1) et A(P2) font 
entre elles un angle a. La vitesse relative Vr, pour une valeur de a, est 
égale à : 

Vr = )/VI 2 + V2= - 2 VI V2 cos a (3). 

Dans un cas général, a peut prendre n'importe quelle valeur comprise 
entre 0 et 2 7r, toutes les directions étant équiprobables. 

La probabilité d’avoir une valeur comprise entre a + e et a — e est : 

da 


2 77 

Le phénomène étant continu (« peut varier de 0 à 2 de façon continue), 
on peut intégrer l’équation 3 pour avoir la valeur de Vr. 

On a : 


1 

Vr =- 

2 7r 


2 — 


Vr da 


r 2 tt 

1 __ 

- 1 VI 2 + V2 2 - 2 VI V2 cos ada 

2 7r 

9J 


0 


0 
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On peut simplifier cette équation en posant : 

2 VI V2 

A =- 

VI 2 + V2 2 


L'équation devient : 

1 _ 

Vr =- 1/ VI 2 + V2 2 

2 7T 


| 1 — A COS a da 


( 4 ). 


r 2 tt 

Résolution de l'intégrale \ (1 —A cos a ) 1/2 da. 

J 0 

Notons d'abord que cette intégrale peut être divisée en deux autres 
intégrales: 

2 TT f IT f 2 7T 

( 1 — A COS a) 1/2 da = \ ( 1 — A cos a) 1/2 da + I ( 1 — A COS a ) 1/2 da 

0 J 0 J, 

Posons : a = 2 w — p. 

Il vient : da = — dp. 

cos a = cos (2 ir — p) = cos p. 

f 2 71- f T 

Donc : 1 ( 1 — A cos p ) 1/2 dp = 2 \ 

JO J 0 

La résolution de cette intégrale est impossible par les moyens classi¬ 
ques. Cependant, elle peut être calculée en considérant qu'elle est égale 
à une somme finie des valeurs particulières de l'intégrale pour des valeurs 
d'a extrêmement petites, variant de 0 à «•. 

Nous avons programmé cette équation en langage machine « Fortran » 
en donnant à a des valeurs successives de 0 à ir, par pas de tt/ 100. Le 
calcul est effectué sur un ordinateur IBM 7040 (3). 

En outre, pour l'utilisation pratique de l’équation (2), il importe de 
connaître Vr pour le plus grand nombre possible de valeur de A. On 
peut facilement en calculer un grand nombre car A est compris entre 
0 et 1. En effet : 

(3) Travail effectué au Centre de Calcul et de Traitement de l'Information de l'Uni¬ 
versité de Liège. 
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A > 0 car les vitesses sont toujours positives ou nulles 

A < 1 car VI 2 4- V2- > 2 VI V2 

en effet : 

VI 2 + V 2 — 2 VI V2 = (VI — V2) 2 > 0 (un carré est toujours posi¬ 
tif ou nul). 

De la même façon que précédemment, nous avons programmé l’équa¬ 
tion (5) en faisant varier A de 0 à 1 par pas de 0.01. 

1 Ç 2 TT 

S =- 1 ( 1 — A COS a) l/ - da (5) 

2 77 J 0 

Les valeurs obtenues sont reproduites dans le tableau 1, L'erreur 
maximum effectuée en utilisant ce mode de calcul peut être évaluée. 
En effet, la courbe représentant l’équation est croissante, car : 


( I 1 — A cos a )-A sin « > 0 si 0 < « < 77 

\ / 2 V" 

Calculons l’intégrale par la méthode des rectangles : 


1 99 T 1 0 77 

Si = — S - | 1 — A cos- 

77 0 = 0 100 100 


1 

< — 


1 100 

<— 2 - 


J 0 


77 0 = 1 100 


1 P* 

| 1 — A cos-= S 2 


100 


S, = 


1 99 

S 


100 0 = 0 


V . 


— A cos 


P * 


100 




1 100 1 / 0 77 


100 0=1 


V . 


— A cos-= So 

100 


L'intégrale vaut donc : 


S = 


S, 4- S 2 


2 


+ A 
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Evaluons l'erreur |A| : 

S,-S, 1 /_. _\ , y 

|A| =-=- 1/ 1 + A - 1/ 1 - A < -= 0,007 

2 200 \ / 100 


L'erreur est maximum quand A est grand (proche de 1). Dans le cas 
considéré, l'erreur est de 0,7 %. 


Applications expérimentales. 

Vérifications des hypothèses. 

Dans toutes les équations développées plus haut, nous avons supposé 
que lès animaux se déplaçaient de façon aléatoire. Dans le cas des Inver¬ 
tébrés, cette hypothèse semble vraisemblable. Berthet (1964) a mesuré 
les déplacements de diverses espèces d'Oribatides marqués au Na-_, Ir 192 
Cl 0 . Les parcours observés rappellent nettement les images obtenues lors 
de l’étude des trajectoires de particules animées d'un mouvement brownien. 

Pour vérifier la concordance des valeurs théoriques déduites de l'équa¬ 
tion (4) et les observations expérimentales, nous avons effectué deux 
types d'expériences. 


Expériences sur Gammarus pulex. 

Dans cette expérience, nous avons étudié un cas particulier de l’équa¬ 
tion (4). Considérons, en effet, que la population PI est constituée par 
des individus immobiles (prédateurs à l’affût, par exemple). On a donc: 
VI =0. 

L’équation est alors simplifiée et on a : 

Vr =- V2 

Pour connaître le nombre de rencontre, il suffit de remplacer les 
valeurs demandées dans l'équation (2). Ces valeurs sont toutes connues. 


Dispositif expérimental. 

Les animaux sont placés, une heure après leur capture, dans une 
enceinte expérimentale de 24 X 30 cm contenant de l'eau sur une hauteur 
d'un cm . L'eau est maintenue à une température constante de 10 "C. Sur 
le fond, se trouvent des cercles noirs représentant les animaux immobiles 
de la populations PI. Ces cercles ont un diamètre d'un cm et sont distri¬ 
bués au hasard, à une concentration de 0,11 points par cm 2 . 


DE POPULATIONS ANIMALES HOMOGENES 

























8 


R. R. TERCAES. — RENCONTRES ENTRE INDIVIDUS 


41 36 


L'enceinte est filmée, sur Gevapan 30, pendant la durée de l'expérience, 
à l'aide d’une caméra Pathé Webo M 16 mm, équipée d’un objectif de 
25 mm de focale (Angénieux 1 : 1.4), à la cadence de 8 images/seconde. 

L'examen du film, à la visionneuse, permet de déterminer la vitesse 
des animaux P2 (Gammarus) et le nombre de rencontres entre les 
animaux Pl et P2. 

Le diamètre de perception des Gammarus est donné par la largeur 
de leur corps, dans le sens du déplacement : 0.3 cm. 

Pour augmenter la précision, nous avons examiné chaque animal sépa¬ 
rément. La vitesse individuelle est ainsi chaque fois calculée à chaque 
déplacement. Ce procédé supprime également l’obligation de fractionner 
le temps expérimental en temps infiniment petits dt dans le cas où une 
partie de la population reste immobile. 


Résultats. 

Les résultats obtenus sont repris dans le tableau 2. On constate une 
excellente relation entre les résultats obtenus par l'observation et ceux 
obtenus par le calcul, en ce qui concerne la moyenne des diverses 
expériences. Cependant, les fluctuations individuelles sont assez impor¬ 
tantes ce qui entraîne une valeur non significative dans le calcul d'un 
test de corrélation. Ces fluctuations sont tout-à-fait normales vu le nom¬ 
bre de mesures effectué. La bonne concordance entre les résultats moyens 
calculés et observés prouvent que les animaux se sont bien dirigés au 
hasard et que la formule proposée est parfaitement valable, en ce qui 
concerne le nombre de rencontres moyen. 


Expériences sur Steganacarus magnus. 

L'expérience suivante est effectuée en utilisant les données de Berthet 
(1964, page 12). En effet, cet auteur a publié une figure représentant les 
déplacements de cinq individus de l'espèce Steganacarus magnus, en 
fonction du temps. De cette figure, on peut déduire la vitesse de dépla¬ 
cement moyenne de chaque individu. On peut donc, à partir de ces 
données, effectuer une expérience identique à celle conduite sur les 
Gammarus. 


Dispositif expérimental. 

Sur la figure représentant le déplacement des Oribatides (popula¬ 
tion P2), nous plaçons une surface transparente contenant des points 
noirs représentant les animaux immobiles Pl, distribués au hasard. Les 
points noirs ont un diamètre d'un cm et sont distribués à la concentration 
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TABLEAU 2 

Expériences sur Gammarus pulex. Nombre de rencontres observés et calculés. 


Animal n° 

Vitesse cm/sec. 

Nb de rencontres 
observé/sec. 

Nb de rencontres 
calculé/sec. 

1 

3.4 

0.0 

0.49 

2 

6.4 

0.8 

0.92 

3 

5.2 

0.8 

0.74 

4 

3.0 

1.0 

0.43 

5 

6.0 

1.4 

0.86 

6 

6.0 

1.0 

0.86 

7 

5.8 

0.8 

0.83 

8 

3.8 

1.0 

0.54 

9 

6.7 

0.0 

0.96 

10 

7.0 

2.0 

1.0 

11 

5.7 

0.7 

0.81 

12 

5.0 

1.1 

0.72 

13 

5.7 

1.0 

0.81 

H 

4.5 

0.0 

0.64 

15 

3.7 

0.4 

0.53 

16 

4.8 

0.5 

0.69 

17 

6.5 

0.5 

0.93 

18 

4.7 

1.0 

0.67 

19 

8.0 

1.3 

1.14 

20 

4.3 

0.0 

0.61 

Moyenne 

5.31 

0.765 

0.759 


de 0.11 individus/cm-. La distance de perception des Oribatides est de 
0.1 cm. 

L'examen des deux figures superposées donnent immédiatement le 
nombre de rencontre total. Dix expériences ont été effectuées sur chaque 
tracé. 
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TABLEAU 3 

Expériences sur Steganacarus magnus. 

Nombre de rencontres observés (moyenne de 10 expériences) et calculés. 


Animal n° 

Vitesse cm/lour 

Nb de rencontres 
observé/jour 
(Moyenne 
de 10 exp.) 

Nb de rencontres 
calculé/jour 

1 

3.90 

0.48 

0.47 

2 

5.02 

0.55 

0.60 

3 

3.90 

0.43 

0.41 

4 

4.22 

0.52 

0.51 

5 

3.01 

0.38 

0.36 

Moyenne 

4.01 

0.472 

0.470 


Résultats. 

Les résultats obtenus sont présentés dans le tableau 3. Comme dans 
le cas des Gammarus, il existe une excellente correspondance entre les 
résultats moyens observés et les valeurs déduites des équations. Les ani¬ 
maux se sont donc bien déplacés au hasard. 


DISCUSSION ET CONCLUSION. 

Plusieurs modèles mathématiques ont déjà été proposé pour exprimer 
les variations de populations au sein d'une communauté animale simpli¬ 
fiée, contenant des prédateurs et des proies ou des parasites et des hôtes, 
La méthode classique, celle de Lotka (1925) et Volterra (1939), est de 
faire intervenir des coefficients de natalité et de mortalité, spécifiques 
aux deux niveaux trophiques de la population étudiée. 

Des modèles analogues ont été présentés par la suite (Bodenheimer 
et al., 1952; Ullyett, 1953) où interviennent également des coefficients 
regroupant un grand nombre de phénomènes biologiques. 

A ces modèles déterministes simples, d’autres auteurs ont tenté de 
substituer des modèles stochastiques (Bartlett, 1957; Bailey, 1964). 
La plupart de ces modèles conduisent à des équations insolubles dont 
les solutions ne peuvent être connues que par approximation (Whittle, 
1952). D'autre part, l’application de ces modèles à des situations réelles. 
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par exemple la croissance d’une population de paramécies sur une culture 
de levure et de bactéries, montre qu'il existe rarement une situation écolo¬ 
gique suffisamment homogène pour être exprimée par ce type de modèle 
mathématique (Bartlett, 1957, 1960 discutant les résultats expérimen¬ 
taux de Gause, 1934). 

Ces diverses considérations nous ont amener à reconsidérer le pro¬ 
blème d'un modèle mathématique de populations en partant de notions 
biologiques plus élémentaires. 

Prenons, par exemple, le « coefficient de natalité » utilisé dans presque 
tous les modèles étudiés. Ce coefficient tente de couvrir un éventail impres¬ 
sionnant de phénomènes. Citons en quelques-uns : existence de mâles e t 
de femelles dans une population donnée à un moment donné, rencontre 
entre les mâles et les femelles, accouplement, gestation, ponte des œufs, 
développement des œufs, développement des larves et passage au stade 
adulte. Ces phénomènes ont des durées différentes et peuvent être affec¬ 
tés par divers facteurs extérieurs : destruction des œufs par les prédateurs, 
mortalité des jeunes pour des causes physiques, etc. Il est donc illusoire 
de vouloir représenter une telle masse de processus biologiques par un 
seul paramètre mathématique, même considéré comme variant de façon 
aléatoire. 

Le travail que nous avons entrepris et dont les grandes lignes ont été 
définies au début de cet article a pour but de substituer à divers coeffi¬ 
cients globaux toute une série de paramètres calqués sur la réalité biolo¬ 
gique. Le détail des paramètres ainsi choisis et la justification de ce choix 
seront exposés dans un travail qui sera publié ultérieurement. 

Une valeur essentielle intervenant dans un grand nombre de situations 
biologiques est le nombre de rencontres entre les constituants d'une popu¬ 
lation; par exemple, le nombre de rencontres prédateurs-proies, mâles- 
femelles, adultes-jeunes, etc. 

Les équations que nous avons mis au point (équations 2, 4 et tableau 
1 ) permettent de calculer le nombre de rencontres entre les individus 
appartenant à n populations animales différentes. Elles tiennent compte 
des vitesses de déplacement respectives de chaque composant des popu¬ 
lations étudiées, de leur distance de perception, de leur zone d'influence 
et de leur densité. 

Tout facteur extérieur influençant l'un ou l'autre de ces paramètres va 
modifier le nombre de rencontres. Ainsi, il est bien évident que tout facteur 
modifiant les vitesses de déplacement de l'une ou l’autre catégorie de 
populations va affecter la vitesse relative Vr. La température, la saison, 
le niveau d'alimentation vont intervenir. Il importera donc d’en tenir 
compte dans tout problème étudiant une population pendant une longue 
période. Un raisonnement analogue s'applique aux dimensions des dia¬ 
mètres de perception et d’influence. 

De même tout facteur modifiant la concentration relative d'une ou de 
plusieurs catégories d'animaux va modifier le nombre de rencontres. 
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Il faudra donc fractionner la durée du phénomène étudié en périodes dt 
convenables. Le bilan total pourra être une intégration si le phénomène 
est continu ou une simple somme si le phénomène est discontinu. 

Une équation analogue mais d'utilisation moins générale a été formulée 
par Skellam (1958) dans un problème différent: le calcul de la densité 
d'une population d'oiseaux en considérant le nombre d'oiseaux rencontrés 
par un observateur se déplaçant dans une forêt. Skellam a effectué la 
résolution de l’intégrale exprimant la vitesse relative en utilisant les tables 
de Legendre (1825), méthode moins générale que celle que nous avons 
utilisée (tableau 1). Dans la suite de notre travail, cette résolution précise 
aura une grande importance. 

Pour terminer, il nous faut encore insister sur le fait que la formulation 
proposée ici n'est valable que si les animaux considérés se déplacent au 
hasard pendant l’espace de temps considéré. Une telle hypothèse de tra¬ 
vail semble raisonnable dans le cas des Invertébrés et est confirmée par 
nos expériences sur Gammarus et Steganacarus. 

Les équations que nous proposons constituent donc un outil de travail 
intéressant dans le cas des études de populations car elles permettent de 
quantifier un grand nombre de phénomènes biologiques. Nous nous ser¬ 
virons de ces équations dans le programme général défini au début de cet 
article. 


RESUME. 

Présentation détaillée d’une équation exprimant le nombre de rencon¬ 
tres (0 12 ) entre les individus de deux populations animales (1 et 2), en 
fonction de divers facteurs : 

d , 2 — n m D Vr 

n = nombre d'individus dans la population 1 
m = nombre d'invidus/cm 2 dans la population 2 
D = somme des diamètres d’influence ou de perception Dl et D2 
Vr = vitesse relative des individus 1 par rapport aux individus 2. 
Extension à un système à n populations. 

Exemples d'expériences. 


SUMMARY. 

Detailed présentation of an équation giving encounter number ( 6 V .) 
between the individuals of two animal populations (1 and 2), in terms 
of several parameters : 

6,2 — n m D Vr 

n = number of individuals in population 1 
m = number of individuals/cm 2 in population 2 
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D = sum of influence or perception diameters Dl and D2 of individuals 
1 and 2 

Vr = relative velocity of individuals 1 to individuals 2. 

Generalization to a n population System. 

Examples of experimer.ts. 

Institut royal des Sciences naturelles de Belgique. 
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